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Wave  spectra  were  measured  In  the  nearshore  region  extending  approximately 
900  m from  the  shore  at  the  Island  of  Sylt  (North  Germany)  In  the  North  Sea. 

Seven  gage  stations  were  In  operation  for  a period  of  two  weeks  In  May,  1976 
for  data  collection.  The  field  results  were  presented  and  compared  wlt^ 
numerical  computations  for  energy  transformation  in  shallow  water,  based  on  a 
method  developed  earlier  for  the  same  project.  The  numerical  computation 
agreed  well  with  the  field  data  In  the  range  of  energy-containing  wave 
components  and  low  frequency  components  but  failed  in  high  frequency  range. 

Based  upon  commonly-used  wave  Instability  criteria,  an  equilibrium 
energy  spectral  density  function  in  water  of  finite  depth  has  been  developed. 

This  function  provides  saturation  conditions  on  spectral  density  for  the 
complete  spectrum.  Comparison  with  field  data  is  encouraging  but,  nevertheless, 
inconclusive. 

Combining  field  evidence  and  numerical  results,  an  assessment  has  been  made 
on  the  wave  climate  of  the  test  area.  It  was  demonstrated  that  the  offshore 
bar  located  in  the  test  area  plays  a double  role  as  energy  dissipator  and  trapper. 
It  was  also  shown  that  waves  from  the  N-W  quadrant  are  more  effective  in 
generating  longshore  current  and  thus  in  longshore  material  transport.  . 
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1.  Introduction 

The  utillsetion  of  wave  energy  apectma  as  a etatlatlcal  aeaeure  to  describe 
the  ever-lrregular  ocean  wave  has  becoae»  «ore  or  less,  a standard  practice. 

Armed  with  this  tool,  oceanographic  scientists  have  made  significant  progress  In 
the  past  t«K>  decades  In  the  developsient  of  wave  forecasting  techniques.  The 
endeavor  has  been,  however,  mainly  concentrated  on  deepwater  situations. 
Consequently,  the  deep  water  wave  climate  can  be  assessed  with  certain  degrees 
of  confidence.  To  establish  the  wave  climate  In  the  coastal  zone  as  Is  often 
necessary  for  scientific  and  engineering  purposes  the  deepwater  wave  environment 
needs  to  be  transformed  Into  shoaling  waters.  It  Is  In  this  area  that  both 
knowledge  and  effort  are  conspicuously  lacking,  although,  odd  as  it  may  seem, 
the  transformation  of  monochromatic  wave  train  Is  a well  studied  topic. 

At  present,  there  exists  a few  attempts  In  the  development  of  theoretical 
concepts  and  numerical  techniques  concerning  spectral  transformation  from  deep 
to  shallow  water.  Longuet-^Hlgglns  (1956)  and  Collins  (1972),  for  Instance, 
have  studied  the  theoretical  formulation  on  wave  spectrum  by  considering  the 
conservation  of  energy  between  adjacent  wave  rays  much  the  same  way  of  the  wave 
refraction  in  monochromatic  wave  trains.  Based  on  Longuet-^ggins ' formulation, 
Karlesson  (1969)  developed  a method  to  compute  spectrum  transformation  over 
parallel  bottom  contours.  Collins  (1972)  extended  Che  work  to  the  inclusion 
of  bottom  friction.  Krasitkly  (1974)  has  derived  explicit  analytical  solutions 
for  the  spectrum  transformation  over  two-dimensional  parallel  bottom  contours 
More  recently,  Shiau  and  Wang  (1976)  developed  a numerical  technique  to  compute 
the  spatial  transformation  of  wave  spectrum  over  irregular  bottom  topographies. 
Their  work  promises  to  be  an  effective  method,  though  it  has  not  been  critically 
tested  for  lack  of  detailed  field  data  as  yet. 
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In  May,  1976,  a brief,  but  intensive,  field  work  was  carried  out  at 
the  Island  of  Sylt  (North  Germany)  in  the  North  Sea  to  collect  information  on 
nearshore  wave  environment  and  beach  process.  One  of  the  objectives  was  to  better 
the  understanding  of  nearshore  wave  energy  transformation  and  to  test  validity 
and  limitation  on  the  numerical  technique  developed  by  Shiau  and  Wang.  The 
material  covered  in  this  report  deals  specifically  with  this  objective. 

The  field  work  covered  a period  of  approximately  two  weeks.  Wave  spectra 
were  measured  at  stations  Installed  in  water  ranging  from  a 7-meter  depth  to 
less  than  1 meter  near  the  shoreline  at  low  tide.  During  this  period,  the  wave 
conditions  were  considered  normal  for  this  area  with  heights  up  to  1.5  to  2 meters. 
The  winds  varied  in  speed  from  as  low  as  3 m/s  to  more  than  12  m/s;  the  wind 
directions  covered  both  onshore  and  offshore. 

These  field  results  were  compared  with  numerical  predictions  using  the 
recorded  offshore  data  as  input.  During  the  course  of  the  study,  the  numerical 
procedure  was  further  improved  to  relax  the  boundary  condition  requirements  and 
to  shorten  the  computational  time. 

An  equilibrium  spectral  density  function  has  been  developed  for  the  condition 
of  water  of  finite  depth.  This  function  provides  information  on  the  saturation 
condition  of  energy  density  for  the  complete  spectrum.  Finally,  combining  field 
evidence  and  numerical  results,  an  attempt  was  made  to  provide  an  assessment  of 
the  wave  climate  of  the  test  area. 

2.  Theoretical  Background 
Transformation  of  Energy  Spectrum 

In  Its  very  general  form,  the  total  energy  of  a random  seaway  can  be 
expressed  as  the  summation  of  energy  content  in  the  wave  number  and  frequency 
space  such  that 


E = i x(^»n)  dk  dn 


where  E is  the  total  energy;  x(k,n)  is  the  spectral  density  function;  k is 
the  wave  number  and  n is  the  wave  frequency. 

Reduced  spectra  can  be  obtained  from  x(k,n)  by  integration  over  n and 
over  k.  Thus,  we  have 


x(^*n) 

n 


= xCE.n) 


where  F(k)  and  $(n)  are  known  as  the  wave  number  and  wave  frequency  spectral  den- 
sity functions,  respectively. 

Since  the  wave  number  spectrum  can  be  interpreted  as  a directional  fre- 
quency spectrum  owing  to  the  dispersion  relation  between  k and  n,  the  following 
relationship  exists  (Phillips,  1966) 


^(n,e)dn  - [k  F(k,0)  dk) 


k » G(n) 


where  k is  the  magnitude  of  k and  6 is  the  direction  of  propagation  and  k » G(n) 
is  the  dispersion  relationship.  Let  subscript  0 denote  the  reference  state 
(deepwater  condition  is  usually  conveniently  selected),  the  transformation 
equation  of  frequency  spectrum  is  established. 


[k  F(k.e)dk),  . 
Uk  F(k.9)dkl, 


(t(n,e)  dn]g 
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For  steady  state  condition,  the  frequency  n is  invariant  to  space,  the  above 
equation  simplifies  to 


♦(n,e) 


tk  F(k,9)dh],^  . 
[koFo(k,6)dk)^  . G^(„) 


(6) 


Since  F(k,d)dk  approximates  the  incremental  of  wave  energy  over  the  wave 
number  band  dk  at  a directional  angle  6,  it  is  plausible  to  write 


F(k,e)dk  - dE(k;e)  (7) 

with  6 as  a parameter  here. 

Substituting  Eq.  (7)  into  Eq.  (6)  results  in 


sfn,9) 


(k  dE(k;8)Jj.  , 

tk^dE^(k;e))k 


G(n) 

■ G„(n) 


♦o (n,0) 


(8) 


Therefore,  to  compute  the  spatial  variations  of  t,  we  must  first  establish  the 
spatial  variations  of  k and  dE(kje). 


Wave  Number  and  Wave  Energy  Transformation 

The  basic  equation  to  compute  the  spatial  variations  of  k is  the  conservation 
equation  for  the  wave  number,  which,  for  steady  cases,  reduces  to  the  simple 
form 


or 


7n  - 0 

n ■ const  ■ o 

o 


In  the  presence  of  a steady  current  u,  we  have 


(9) 

(9-a) 


[-  For  regions  where  both  d and  u very  slowly  over  distances  of  the  order  of  a tmve- 

^ length,  the  following  dispersion  relationship  exists, 

r 


Coablnlng  Eqs. 


2 

a " gk  tanh  kd 
(9),  (10)  and  (11)  results  In 
1/2 

(gk  tanh  kd)  u k cos6  + v k sin  6 


(11) 


(12) 


for  the  coordinate  system  shown  in  Fig.  1.  Here  u and  v are  current  velocity 
components  In  the  x and  y directions  respectively;  6 is  the  wave  approaching 
angle  measured  from  the  positive  x^axls. 

Since  It  is  known  that  the  wave^number  vector  In  space  Is  irrotatlonal. 
It  follows 


I 


3(k  cos  9)  _ 3(k  sin  0)  _ j 

3y  ” 3x  ''  ' j 

For  given  spatial  distribution  of  u,  Eqs.  (12)  and  (13)  enable  us  to  compute  '• 

k and  9 In  the  wave  region  provided  the  boundary  condition  Is  specified.  * 

To  determine  the  variations  of  dE(k;6),  the  assumption  is  made  that  the  ] 

wave  energy  associated  with  a certain  frequency  band  stays  within  this  band  so  1 

that  the  principal  of  energy  conservation  is  applicable  to  each  individual 

frequency  band  under  consideration.  This  assumption  denies  energy  transfer 

among  different  frequencies  and  Is  compatible  only  when  nonlinear  effects  are 

negligible.  Following  the  energy  equation  presented  by  Longuet-Hlgglns  and 


) 
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Stewart  (1960,  1961)  for  the  fluctuating  motion  of  waves  with  a superimposed 
current  system,  we  have 


[dE(k)(u  + Cg  cos  e)]  + ^ (dE(k)(v  + Cg  sin  6)]  + 

+ S + S - -E 

xy  9x  yy  9y  yx  3y 


where  Cg  is  the  group  velocity  expressed  as 


„ - ^ - 1 n + 2kd 

g C 2 * sinh(2kd)^ 


C is  the  wave  celerity 

S ^ are  the  radiation  stresses  such  that 
aB 

^xx  “ dE(k)[(2ng  - |)  cos^e  + (n^  - |)  sin^0]  » dE(k) 

S “ dE(k)[(2  n - 4)  sin^0  + (n  - 4)  cos^6]  - dE(k)  a 
yy  o ^ s ^ yy 


S - S » n„  sin  (26)  - dE(k)  o^ 

xy  yx  2 g ' ' ' ' xy 


£ is  the  rate  of  energy  dissipation  per  unit  area* 


With  known  9,  Eq.  (14)  permits  the  computation  of  the  -spatial  variations  of  dk(E) 
for  given  boundary  conditions  and  the  dissipation  function.  For  the  present 


computation,  1 is  taken  to  be  equal  to  zero. 


3.  Numerical  Procedures 


Based  upon  Eqs.  (12),  (13)  and  (14),  Noda,  et  al.  (1974)  developed  numerical 
procedures  to  compute  the  wave  number  and  wave  height  as  functions  of  space 
variables  for  beaches  with  periodical  bottom  variations  in  the  longshore  direc- 
tion. Shiau  and  Wang  (1976)  extended  the  procedure  to  compute  the  spectrum 
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transformation  with  the  boundary  conditions  somewhat  relaxed  to  accow>date 
irregular  but  slowly  varying  bottom  profiles.  In  the  present  computation, 
their  method  is  further  improved  to  relax  the  requirements  of  boundary  conditions 
and  to  Improve  the  computational  efficiency.  The  procedures  are  briefly  des- 
cribed here. 

To  facilitate  numerical  computation,  the  governing  equations  are,  out  of 
necessity,  to  be  transformed  into  finite  difference  form.  A mixed  forward, 
backward  and  central  differences  scheme  are  used  to  minimize  the  boundary 
restrictions  at  the  expense  of  nonuniformity  of  the  order  of  errors.  With 
reference  to  the  grid  schemes  shown  in  Pigs.  2 and  3 the  key  equations  are 
suomaarlzed  here: 

(1)  The  wave  number  and  wave  angle  are  solved  through  numerical  iterations 
of  the  vave-nuaber  conservation  equation  and  the  irrotational  wave-number 
equation.  The  conservation  equation  [Eq.  (12)], 


[g  k,j  tanh  (k^j  d^j)]  ' + k^j  coa  k^^  sin  6^^  - 


la  solved  through  a Newtonian  Iterative  technique  defined  as: 


V V .1!^ 

new  “ old  f ' 

new 


where  f*(k)  is  the  first  derivative  of  f(k)  which  is  defined: 


f(k)  ■ gk  tanh  kh  - - uk  cos  e - vk  sin  6] 


The  iteration  is  considered  satisfactory  when 


Ik  - k . .j  < 0.001  |k^ 
' new  old’  — ' n 
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t 


The  computation  is  carried  to  the  extent 

le  - 6 , < 0.001 

' new  old'  — ' new' 

Since  forward  difference  is  used  in  Che  x-dlrectlon  and  backward  difference  is 
used  in  the  y direction,  6 needs  to  be  specified  at  the  offshore  boundary  and 
at  the  x-axis.  The  offshore  boundary  condition  is  the  input.  The  condition  of 
the  x-axis  could  be  specified  in  such  a manner  that  0^^^  - 9^^  without  introducing 
significant  error  if  the  longshore  bottom  variation  in  the  vicinity  of  the 
X'axls  is  slow. 

The  above  numerical  scheme  becomes  unstable  when  0.  The  exact  stability 

criterion  cannot  be  easily  determined.  However,  using  the  null-current  condition 
as  a guideline,  one  arrives  at  the  following  condition: 


> Tan 
^X 


where  9 is  the  maximum  expected  wave  angle.  For  example,  if  the  maximum  wave 

IB 

angle  is  expected  to  be  leas  than  60”,  the  ratio  of  dy/dx  should  be  larger  than 
1.7.  The  velocity  and  depth  gradients  required  in  the  computation  are  determined 
by  central  difference,  i.e.. 


IJ  dx, 


, Ai_  . ^i+l.J  ~ 


where  ^ represents  u,  v or  d and  Xj^  is  either  x or  y. 

The  gradients  at  the  boundary  are  taken  as  equal  to  xero  at  offshore  and  both 
sides  and  as  constant  at  the  shoreline,  that  is 


These  boundary  conditions  can,  of  course,  be  easily  altered  to  fit  Individual 


cases. 
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(2)  The  wave  energy  equation  In  finite  difference  form  la, 

[dE(k)]^^^^j^(v  +Cg  sine)  Ax  - [dE(k)  (u  + CgCOse)Ay 


ij 


(v  + Cg  slne)Ax  ” CO®®)  Qj^j  AxAy 


(19) 


where  “ W ^ cose)  ^ ®ln®)  + 


_ _ _ j.  „ /®v  . 3uv  . ^ 3v 

°xx  ax  ®xy  ^ax  ay^  ®yy  3y 


are  as  defined  In  Eq.  (16). 


Here  all  the  values  of  u,  v and  6 are  evaluated  at  (1,J)« 
Eq.  (19)  Is  solved  through  a row  by  row  relaxation  until 


H - H , , 
' new  old 


< 0.01  |h  I 

— ' new' 


(20) 


Like  the  wave  angle  solution,  the  finite  difference  scheme  Is  forward  in  the 
x-dlrectlon  and  backward  In  the  y-directlon.  The  boundary  values  are  also 
defined  in  a similar  manner  In  that  the  shoreward  boundary  Is  specified  as 
input  and  the  values  on  the  x-axls  Is  taken  the  same  as  the  next  column  (E. 


11 


E,2). 


On  the  shoreward  direction,  the  wave  will  eventually  break,  the  dissipation 
term  e In  Eqs.  (14)  can  no  longer  be  neglected.  If  dissipation  Is  not  allowed 
for,  the  wave  energy  density  will  keep  Increasing  and  the  numerical  solution  will 
eventually  become  unstable.  Therefore,  a constraint  must  be  built  In  to  restraint 
super  saturation  of  energy  flux.  Since  the  real  mechanism  of  energy  dissipation 
is  not  well  understood.  An  esiplrical  criterion  based  upon  total  wave  energy  In 
terms  of  significant  wave  height  is  employed  here,  such  that  the  height  Is 
restricted  by  the  following  condition  (Dlvoky,  et  al.,  1970) 


(—■)  - 0.12  t.nh  2»  (^) 

’•S  b *-S  b 


(21) 
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where  the  subscript  b refers  to  breaking  condition  and  Hg  Is  related  to  the 
total  wave  energy  in  the  conventional  expression: 

Hg  - 2.828  (22) 

This  criterion  la  to  be  further  examined  in  later  section  when  we  have  the 
field  data  in  hand. 

The  Input  requirements  Includes  a grid  of  water  depth  augmented  by  the 
tidal  condition.  The  grid  size  should  be  selected  In  accordance  with  the 
criterion  set  In  Eq.  (18).  Strictly  speaking,  directional  spectra  of  the 
form  ^g(nt@)  should  be  prescribed  along  the  seaward  boundary.  Since  our  present 
knowledge  on  directional  spectrum  la  quite  Inadequate,  the  program  will  also 
accept  one-dlmenslonal  spectra  as  the  Input.  The  directional  Information  should 
then  be  provided  as  predominant  wave  angles  at  each  grid  point  along  the  outer 
boundary. 

4.  Field  Program 

The  field  experiments  were  conducted  at  the  Island  of  Sylt,  West  Germany. 

The  site  location  and  area  hydrograph  are  shown  in  Fig.  4.  The  site  Is  exposed 
to  the  North  Sea  and  is  characterized  by  the  existence  of  anunderwater  transverse 
ridge  (a  pronounced  sand  bar)  running  parallel  to  the  shore  approximately  500  m 
from  the  shoreline.  This  underwater  sand  bar  is  not  stationary  but  fre- 
quently Is  shifting  its  location  from  season  to  season  and  sometimes  after  a major 
storm.  Apparently  Influenced  by  this  ridge,  the  nearshore  flow-pattern  and  the 
resulting  nearshore  bottom  topography  become  rat^  complicated.  Shallow  scouring 
pits,  of  ten  appearing  In  eddy-llke  pairs,  dotted  the  area  at  uneven  spaces.  Since 
the  normal  tidal  range  Is  2 m In  this  area  and  the  normal  water  depth  over 
the  ridge  Is  only  of  the  order  of  3 m,  waves  often  break  over  the  ridge  and  reform 
afterwards. 


To  measure  the  change  of  wave  characteristics  over  this  region,  seven  wave 
gages  were  Installed  at  locations  as  shown  in  Fig.  3.  Gages  1 and  2,  which 
were  bottom-oaounted  in  relatively  deep  water,  were  echo-transceivers  as  developed 
by  Fuhrentholz  Laboratory  (See  Reference).  These  echo  sounders  were  capable  of 
measuring  water  surface  variations  to  + 0.5  cm  for  water  depth  up  to  90  m.  The 
shallow  water  gages  (3  to  7)  were  staff-mount  pressure  transducers  of  Type  MDS76 

2 

as  manufactured  by  H.  Malhak  AG,  Hamburg  with  a pressure  range  of  0-1  kg/cm  . 

Gages  1 to  5 were  Installed  perpendicular  to  the  shoreline  to  measure  the  wave 
spectra  transformation  over  the  ridge  and  in  shoaling  water.  Gages  5 to  7 were 
mounted  parallel  to  the  shoreline  to  measure  mainly  the  wave  approaching  angle. 

Wave  data  were  collected  from  Ilay  1 to  May  15,  1976  at  selected  times  to 
cover  various  combinations  of  tidal  stages,  waves  approaching  angles  and  wind 
conditions.  Data  were  transmitted  to  shore  station  by  underwater  cable  and  were 
recorded  simultaneously  on  FM  analog  magnetic  tape  and  strip  charts.  During  the 
course  of  wave  measurement,  currents  were  also  measured  with  two-component  electro- 
magnetic current  meters  mounted  at  Gage  stations  2 and  5.  Variable-depth 
drogues  were  also  deployed  at  times  to  determine  the  current  condition.  Descrip- 
tion of  instrumentation  and  recording  devices  was  presented  by  Fuhrbbter  and 
Siisching  (1974) . 

3.  Data  Analysis 

Wave  energy  spectra — Spectral  analysis  was  performed  using  FFT  (Fast  Fourier 
Transform)  technique  developed  by  Cooley  and  Tukey  (1963)  with  a Hanning  Window 
(Black  and  Tukey,  1958)  for  data  smoothing.  All  the  data  sets  were  sampled 
uniformly  at  a time  Interval  At  • 0.25  seconds  to  a total  data  points  N * 2048 


(or  512  seconds  data  period).  These  values  were  chosen  to  cover  adequately  waves  with 


Gage  locations  and  nearshore  hydrograph 
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periods  ranging  from  two  to  nine  seconds  which  were  determined  (as  the  range  of 
energy  containing  tmves)  through  visual  examination  of  the  data  record  and  the 
irlnd  condition  during  measurements.  A section  of  data  sample  is  shown  in  Fig.  6. 


CHART  SPEED  • 2 


After  several  trial  runs,  a degree  of  freedom  DF  « 40  was  selected  as  the  optimum 
to  yield  sufficiently  smooth  spectral  curves  to  facilitate  comparisons  with 
theoretical  prediction,  yet,  at  the  same  time  to  preserve  sufficient  detail  to 
truthfully  portray  the  actual  conditions.  Figure  7 shows  a comparison  between  the 
perlodogram  and  the  smoothed  spectrum. 

Data  collected  In  shallow  water  by  gages  3 to  7 were  actually  wave  pressure 
Instead  of  the  true  water  surface  fluctuations.  To  derive  the  energy  spectra 
from  the  measured  pressure  signals,  the  depth  attenuation  effect  needs  to 
be  compensated  according  to  the  relation: 


S(f) 


. cosh  kd 
^cosh  k(d-t) 


Sp(0 


(23) 


where  S^Cf)  is  the  pressure  spectral  density  as  obtained*  ^ Is  the  depth  of 
the  transducer  measured  from  me'^n  water  level. 


1 


Fig.  7:  Comparison  between  periodogran  and  smoothed  spectrum 


This  correction  factor  is  based  upon  linear  wave  theory  and  it  tends  to  over 
compensate  high  frequency  components  which  were  more  or  less  contaminated  by 
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noise.  An  attempt  was  made  to  apply  this  correction  factor  to  data  analysis  and 
was  later  abandoned  because  such  corrections  did  little  to  wave  components  in  the 
energy  containing  range  but  grossly  over  anplifying  the  high  frequency  components. 

Wave  approaching  angle — To  calculate  the  transformation  of  wave  spectra, 
one  requires  to  prescribe  the  one-diraenslonal  wave  spectra  in  frequency  or  period 
domain  as  well  as  the  directional  characteristics  at  the  outer  boundary.  Based 
upon  the  computational  method  suggested  here,  the  directional  Information  should 
be  expressed  in  terms  of  dominant  angle  of  approach  and  if  desirable  a predesig- 
nated spreading  function.  Since  the  present  field  set-up  yielded  no  directional 
information  at  the  deep  water  sites,  an  alternative  method  was  used  which 
utilizes  the  Information  obtained  from  the  shallow  water  gage  arrays.  As  illus- 
trated in  Fig.  8 when  a wave  front  passing  by  a cluster  of  three  wave  gages  set 
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apart  with  known,  but  small  distances,  the  following  relations  are  approximately 
valid : 


C 

C 


12 

23 


"12 


'■23 


12 


23 


cos  6 
sin  6 


(24) 


where  C Is  the  mean  wave  celerity;  S Is  the  horizontal  distance;  At  Is  the 
time  span  for  wave  front  travelling  from  one  wave  staff  to  the  other  and 
6 Is  the  wave  approaching  angle.  The  Indices  1,  2 and  3 Indicate  the  gage 
numbers.  The  above  equation  when  solved  for  angle  6,  yields 

e - Tan‘^(|i^)(^)(-^)]  (25) 

®23  'l2 

For  shallow  water  condition,  we  further  assume 


where  d and  H are  the  mean  water  depth  and  wave  height  respectively.  Sub- 
stituting Eq.  (26)  Into  Eq.  (25)  one  obtains 


(27) 

d,2  H,,  ^23  ^42 


Since  d and  S are  known,  H and  t can  be  measured  off  directly  from  the  wave 
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records  as  Illustrated,  the  angle  of  approach  6 can  be  computed.  The  deepwater 

wave  angle  9 was  then  determined  by  conventional  wave  refraction  coflq>utation 
o 

for  the  dominant  wave  component  or  components. 

6.  Field  Results 

A selected  set  of  data  were  presented  here  to  demonstrate  the  spectral 
transformation  under  various  combinations  of  wave  Incident  angle,  wind  condition 
and  tidal  stage.  Table  1 summarizes  the  cases  reported.  The  wind  record  during 
the  measuring  period  is  displayed  in  Table  2.  Figures  9 to  15  plot  the  wave 
spectra  as  measured  at  various  gage  stations  along  the  control  line  for 
each  case  reported  herein.  The  estimated  wave  and  swell  approaching  angles 
and  the  wind  directions  are  shown  in  Fig.  16,  The  wind  direction  is  defined 
in  terms  of  the  dominant  direction  over  an  8<>hour  period  prior  and  up  to  the 
measuring  time.  The  directional  variations  over  this  period  are  also  included  in 
Fig.  16.  Of  the  seven  cases,  1,6  and  7 are  mainly  wind  waves  as  clearly  revealed 

by  the  shape  of  the  spectra;  the  wave  approaching  angles  all  coincide  with 

the  wind  direction.  Cases  3 and  5 are  dominated  by  swell;  the  wind  direction 
and  wave  angle  bears  no  relation.  Finally,  cases  2 and  4 represent  mixed 
swell  and  sea  condition  such  that  the  wave  angles  are  also  mixed.  In  Case  2, 

for  instance,  the  directions  of  swell  and  sea  are  alntost  perpendicular  to 

each  other. 

The  offshore  shoal  apparently  plays  an  Important  role  In  the  wave  energy 
transformation  mainly  because  the  water  depth  is  only  2.5  m at  the  crest  of 
the  shoal  (see  Fig.  S ) but  the  tidal  range  is  of  the  order  of  + 1 m.  Thus, 
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during  low  tide,  it  acts  like  a submerged  breakwater.  As  the  waves  break  over 
the  shoal,  a significant  portion  of  wave  energy  is  being  dissipated.  Wave  com- 
ponents in  the  energy-containing  range  take  the  heaviest  toll  while  the  long  wave 
components  are  less  affected.  The  high  frequency  components  remain  largely 
intact  and  may  actually  show  some  gain,  probably  through  reclaiming  partially 


the  energy  lost  by  the  main  wave  components.  There  is  no  indication  that  waves 
will  reform  after  breaking.  The  dominant  wave  components  will  not  recover  their 
energy  and  their  form.  The  appearance  of  the  waves  as  well  as  the  energy 
distribution  among  various  frequencies  are  quite  dissimilar  on  the  two  sides  of 
the  shoal.  For  nonbreaking  situations,  the  shoal  tends  to  attenuate  the  peaks 
of  the  spectra  for  onshore  wind  conditions.  When  the  wind  is  offshore,  on  the 
other  hand,  the  peak  is  significantly  amplified  (Case  5)  due  to  the  existence 
of  the  shoal.  This  is  probably  due  to  the  combined  effects  of  opposing  current 
and  the  fact  that  energy  generated  by  winds  inshore  of  the  shoal  is  being 
entrapped.  This  case  will  be  further  analyzed  later. 

In  general,  the  energy  spectra  will  peak  further  when  approaching  the 
shore  before  finally  breaking.  The  total  energy,  however,  tends  to  diminish 
somewhat  in  comparison  with  the  offshore  waves.  Other  measurements  under 
storm  conditions  (Biisching,  1975  ) in  the  same  area  revealed  similar  results. 

For  the  wind  wave  conditions,  the  JONSWAP  spectra  (Hasselmann,  et  al., 

1973)  which  assume  the  following  form  quite  adequately  describe  the  offshore 
wave  conditions: 

-(f-fj^ 

E(f)  - ag^(27T)"^  exp(-  I (|-)*^)  ^ 2 (28) 


0.007  for  f < f 
— m 


0.009  for  f < f 
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Containing  five  free  parameters  f^,  a,  y»  and 

Here  f^  represents  the  frequency  at  the  maximum  of  the  spectrum  and  the  parameter 
a corresponds  to  the  usual  Phillips  constant.  The  remaining  three  parameters 
define  the  shape  of  the  spectrum:  y Is  the  ratio  of  the  maximum  spectral 

energy  to  the  maximum  of  the  corresponding  Plerson-Moskowltz  (1964)  spectrum 

E^(f)  - og^  (2it)"*  exp(-  I (|-)'*]  (29) 

m 

with  the  same  values  of  a and  f : a and  a,  define  the  left  and  right  side 
, m a b 

widths,  respectively,  of  the  spectral  peak. 

Although  the  JONSWAP  spectra  has  5 free  parameters,  among  them,  the  peak 
enhancement  factor  y is  the  only  critical  one  that  makes  It  differ  from  the  P-M 
spectrum.  Figures  17  to  20  offer  some  comparison  among  the  data,  the 
JONSWAP  spectra  and  the  P-M  spectra. 

7.  Numerical  Results  and  Comparison 

One  of  the  main  purposes  of  this  study  is  to  compare  the  field  results  with 
the  predictions  based  on  the  numerical  technique  developed.  Before  the 
computer  program  is  applied  to  the  actual  field  condition,  a number  of  simple 
cases  were  computed  to  Illustrate  the  general  Influence  of  bottom  topography 
on  spectrum  transformation.  First  of  all,  a case  against  rhythmic  hydrograph 
and  shoreline  (Fig.  21  ) Is  Illustrated  In  Fig.  22  . The  results  clearly 
Indicate,  as  expected,  that  the  energy  converges  to  the  convex  section  of  the 
shoreline  (Point  A)  and  diverges  from  the  concave  section  of  the  beach.  A 
case  is  then  shown  with  a two-dimensional  hydrograph  that  consists  of  a 
section  of  negative  slope,  l.e.,  the  equivalent  to  an  offshore  bar.  The  result- 
ing spectrum  transformation  Is  shown  In  Fig.  23«  For  this  case,  the  energy 
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17:  Comparison  among  field  data,  JONSWAP  and  P-M  spectra  (Data 

from  May  3,  1976  at  15:45) 
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19:  Comparison  among  field  data  JONSWAP  and  P-M  spectra  (Data 

from  May  11,  1976  at  10:15) 
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spectrum  first  peaks  up  on  the  positive  slope  before  reaching  Che  bar;  the 
peak  then  diminlahes  in  magnitude  due  to  the  effect  of  the  negative  slope  before 
it  increases  again  when  the  %raves  proceed  towards  the  shore.  Finally*  (beyond 
the  breaking  point)  the  wave  breaks  and  the  spectrum  diminishes  both  in  inten- 
sity and  in  total  energy. 

Computations  for  the  actual  field  condition  were  carried  out  for  tvo 
different  categories.  In  the  first  category*  wave  spectra  as  measured  at  S2 
(225N  from  the  baseline)  were  used  as  inputs  to  compute  their  transformations 
into  nearshore  zone  close  to  the  bench  (S3  to  S3).  The  bottom  hydrograph  used 
in  this  set  of  computations  was  based  upon  soundings  taken  during  the  period 
of  field  measurement.  The  input  bottom-depth  chart  is  shown  in  Fig,  24,  Here 
the  grid  sizes  are  Ax  • 5m  and  Ay  ■ 20m;  the  corresponding  grid  points  are  36 
and  11,  respectively.  Since  the  input  wave  spectrum  was  available  at  only  one 
grid  point  where  the  spectrum  was  actxially  measured,  the  input  along  the  rest 
of  the  offshore  boundary  need  to  be  estimated.  This  was  accomplished  by  first 
carrying  the  known  spectrum  as  measured  back  to  the  deepwater  condition.  This 
deepwater  condition  is  then  assumed  to  be  homogeneous  over  a horizontal  dimen- 
sion large  enough  to  cover  the  width  of  the  offshore  boundary  and  thus  allows 
to  be  used  as  the  input  to  compute  the  boundary  condition  at  the  other  grid 
points.  This  scheme,  of  course,  allows  only  one  exact  match  at  the  boundary 
point  where  the  spectrum  was  actually  measured. 

Results  of  three  representative  wave-swell  conditions  are  presented  here; 
they  are  swell  dominated  condition  (Case  3),  swell-wave  combined  condition 
(Case  4),  and  wind-wave  dominated  condition  (Case  6),  Comparisons  between  the 
measured  and  computed  spectra  are  summarized  in  Figs.  25  to  27,  respectively, 
representing  Cases  3,  4 and  6.  In  general,  the  prediction  appears  quite  reason- 
able for  wave  components  in  the  energy-containing  range  and  low  frequency  range; 
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In  the  high  frequency  range,  the  prediction  becones  less  desirable.  Based  on 
this  observation,  the  wave  conponents  In  the  energy~contalnlng  range  apparently 
are  the  aost  stable  ones  and  are  affected  aalnly  by  shoaling  and  refraction 
In  their  transfomatlon  prior  to  breaking.  The  slight  over  prediction  at  the 
spectral  peak  Is  Dost  likely  due  to  the  fact  that  the  nearshore  wave  data  Is 
actually  pressure  rather  than  surface  variations  and  that  no  correction  has 
been  made  to  compensate  this  attenuation.  For  Case  3,  the  prediction  was 
actually  carried  Into  the  surf  zone.  The  predicted  spectrum  after  breaking 
appears  quite  reasonable  (Fig.  25)  both  In  order  of  magnitude  and  In  shape. 

This  limited  result  seems  to  support  the  breaking  criterion  expressed  In  Eq.  (21). 

A number  of  factors  could  contribute  to  the  Irregularity  In  the  high 
frequency  range.  One  possible  source  of  error  Is  that.  In  this  range,  the 
variations  In  spectral  density  as  well  as  directions  of  propagation  could  be 
considerable  and  that  results  from  one  point  measurement  might  be  Insufficient 
to  provide  a reasonable  estimation  of  the  Input  along  the  complete  outer 
boundary.  Also,  the  local  generation  and  dissipation,  which  are  most  likely 
to  affect  wave  components  In  this  range,  are  not  Included  In  the  model.  There 
are  at  least  some  Indications  that  such  mechanisms  might  provide  partial 
explanation  on  the  discrepancy  between  measurement  and  prediction.  In  Case  6 
(Fig.  27)  wind  Is  in  the  same  direction  as  wave  propagation,  consequently,  local 
generation  tends  to  over-compensate  dissipation,  thus,  resulting  in  a predicted 
value  smaller  than  measured.  On  the  other  hand,  in  Case  4 (Fig.  26)  when  the 
wind  direction  was  at  a large  angle  with  the  wave  direction,  the  prediction 
becomes  larger  chan  measured.  It  Is  possible  that  in  this  case,  the  locally 
generated  wave  was  running  against  the  predominant  swell  and  thus  resulted  in 
considerable  dissipation  at  high  frequency  range.  A further  examination  on 
the  energy  dissipation  criterion  Is  presented  In  the  subsequent  section. 
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Numerical  computations  were  also  made  for  another  category  using 
measured  data  at  SI  (940M  from  the  baseline)  as  input  to  calculate  the 
corresponding  spectra  at  S2  to  examine  the  wave  transformation  over  the  bar 
Chat  lies  between  SI  and  S2  as  shown  in  Fig.  3.  A much  coarser  grid  system 
was  used  with  x ■ 25  m and  y • 75  m;  the  corresponding  grid  points  are  19 
and  lit  respectively*  Since  no  soundings  were  performed  for  such  a large 
area  (during  the  measuring  period)  the  hydrograph  is  based  upon  the  bottom 
depth  chart  of  1972*  as  supplied  by  the  Leichtwelss-Instltut  fur  Wasserbau, 
Braunschweig,  Germany.  Results  are  presented  in  Figs.  28  and  29.  Figure  28 
represents  a case  where  waves  approaching  from  N-W  quadrant  whereas  the  case 
shown  in  Fig.  29  is  for  waves  coming  from  S-W  direction.  For  these  condi- 
tions, the  prediction  fares  far  better  in  the  former  case  than  the  latter. 

A closer  examination  of  Che  computer  revealed  that  the  spatial  variations 
of  wave  spectrum  are  quite  pronounced  from  grid  to  grid  in  both  longshore 
and  onshore-offshore  directions  for  cases  where  waves  come  from  the  S-W 
whereas  the  contrary  is  true  for  waves  from  the  N-W.  This,  at  least  partially 
explains  why  the  prediction  fares  poorly  for  S-W  waves.  As  a further 
Illustration,  Che  predicted  spectral  density  function  one  grid  distance 
from  the  actual  location  was  also  plotted  in  Fig.  29  for  the  S-W  waves. 

The  spatial  sensitivity  is  apparent  and  is  believed  to  be  the  main 
contributing  factor  to  the  poor  correlation  in  this  case.  The  fact  that 
Che  predicted  values  are  generally  considerably  higher  than  the  measured  for 
S-W  waves  seems  to  suggest  chat  the  offshore  bar  is  a stronger  wave  reflector 
for  waves  approaching  from  this  quadrant  than  that  of  the  N-W  quadrant. 


29:  Comparison  of  measured  and  computed  deep  water  spectra 

(Data  from  May  II,  1976  at  10:15) 
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6«  Equilibrium  Range  in  Water  of  Finite  Depth 

It  la  commonly  accepted  that  gravity  waves  have  a limiting  energy  capacity. 
When  energy  is  added  to  a wave  of  constant  wave  number,  the  wave  height  will 
keep  increasing  until  the  limiting  capacity  is  reached.  Beyond  this  point,  the 
wave  %rlll  break.  In  monochromatic  waves,  this  limiting  condition  is  often  transla- 
kinematlc  terms  such  as  the  limiting  wave  steepness,  the  limiting  height  to 
water  depth  ratio,  or  the  more  fundamental  quantities  of  limiting  water  particle 
acceleration  and  crest  angle.  In  the  case  of  wind  generated  waves,  the  concept 
of  equilibrium  range  is  widely  acclaimed.  In  deepwater,  it  has  been  shown 
through  dimensional  reasoning  ( Phillips,  1958  ) that  in  equilibrium 

range  where  energy  has  reached  saturation,  the  spectral  density  function  takes 
the  following  form: 

♦d(n)  - Bg^  n*^  (30) 

where  8, sometimes  known  as  Philips  constant,  is  determined  from  field  data  as 

-2  -2  -2 
ranging  from  0.8  x 10  to  2.0  x 10  with  an  average  about  1.17  x 10 

In  the  computation  of  wave  energy  spectral  transformation,  one  desires  to 
extend  this  equilibrium  range  into  water  of  finite  depth.  This  equilibrium 
critiera  provides  a physical  limitation  beyond  which  energy  cannot  be 
retained.  The  breaking  criteria  used  in  the  present  model  [Eq.  (21)]  is  based 
on  an  empirical  relation  expressed  in  terms  of  wave  height,  or  equivalently  the 
total  wave  energy.  No  information  can  be  derived  from  this  criterion  as  to  the 
energy  eaturatlon  Holt  in  the  wave  frequency  or  wave  number  domain.  This 
problem  of  energy  saturation  in  water  of  finite  depth  Is  further  examined  here. 


It  has  been  pointed  out  by  Philips  (1966)  that  in  the  context  of  wave  number 
spectra,  the  equilibrium  criterion  expressed  in  Eq.  (30)  Is  consistent  with  the 
occurrence  of  sharp  crests  among  the  waves.  The  limit  sharp  crest  wave  occurs 


when  the  local  fluid  acceleration  at  the  surface  exceeds  a certain  fraction  of 
the  gravitational  acceleration,  g.  If  we  extend  this  physical  reasoning  to 
regions  of  shallow  water  by  assuming  that  a certain  wave  component  of  wave  number 
will  become  saturated  when  the  maximum  vertical  fluid  acceleration  at  the  sur- 
face reaches  a critical  value  say,  ag  with  a a constant  of  proportionality,  the 
saturation  criterion  of  this  wave  number  component  can  then  be  expressed  as 


^ - A^(k^)  canh  k^d  - a (31) 

based  upon  linear  wave.  Here  a^  is  the  vertical  particle  acceleration  and 
A^.(k^)  is  the  limiting  wave  amplitude  associated  with  the  wave  number,  k^. 
Since  A is  related  to  the  wave  energy,  we  further  assume  that 


2 r* 
o L 


g^(k)  dk  - g,(k)  dk 


vhere  g^(k)  is  the  equilibrium  one-dimensional  wave  spectrum  integrated  over 
the  angle  of  spreading  and  y 1*  < constant  of  proportionality.  The  above 
assumption  has  the  physical  meaning  that  the  contribution  to  wave  steepness 

* certain  wave  number  cosq>onent  is  mainly  due  to  waves  of  higher  wave  numbers. 
The  validity  of  this  assumption  may  be  illustrated  in  a much  simplified  situa- 
tion with,  say,  two  wave  components  of  distinct  wave  numbers,  k^  and  k2»  with, 
say,  k2  significantly  larger  than  k^,  superimposing  upon  each  other.  One  then 
expects  that  the  effect  of  the  longer  wave  on  the  shorter  wave  is  mainly  an 
apparent  change  of  mean  water  level  upon  which  the  shorter  wave  rides. 

Neglecting  the  Interaction  between  these  two  wave  components,  the  steepness 
of  the  shorter  wave  can  be  considered  to  be  unaffected  by  the  longer  one. 

The  reverse,  ho%rever,  can  not  be  true  as  any  local  irregularity  as  a result 
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of  the  superimposed  shorter  wave  contributes  to  the  change  of  wave  steepness 
to  the  longer  one.  This  argument  when  extended  to  Irregular  waves  leads 
to  the  condition  expressed  in  Eq,  (32), 

Now  substituting  Eq.  (32)  into  Eq.  (31),  reveals  that  g^(k)  should 
assume  the  follo%rt.ng  form 


2 

- _ ../a.2  kd  sech  kd  + tanh  kd 

®a^“^  (k  tanh  kd)3 


(33) 


This  is  the  functional  relationship  of  the  spectral  density  in  the  equilibrium 

range  k > k,  in  water  of  finite  depth.  Since  2(-)^  is  a constant,  a normalized 
o Y 

function  Qa(k)  can  be  defined  such  chat 


Ga(k) 


^a^^^  _ kd  sech^kd  + tanh  kd 

2(^)^  k^  tanh^  kd 


(34) 


The  function  g^(k)  or  G^(k)  is  the  counterpart  of  Eq.  (30)  for  water  of  finite 
depth.  The  value  of  k g^(k)  approaches  a constant  when  d -*■  « ag  tt  should 
(see  Phillips,  1966).  For  water  of  finite  depth,  however,  k^g^(k)  la  a mono- 
tonlcally  decreasing  function  of  increasing  kd.  This  relationship  in  terms 
2 

of  k G^(k)  vs.  kd  is  shown  in  Fig.  30.  The  variation  of  G^(k)  with  kd  is 
also  computed  and  Is  shown  in  Fig.  31. 

The  constant  In  Eq.  (33)  can  be  estimated  based  upon  the  fact  that  when 
d -»■  « the  functional  relationship  expressed  In  Eq.  (33)  should  approach  Eq.  (30), 


i.e.. 


11m 
d • 


>0 


g^(k)dk  - 


(35) 


where  » »^gk^  is  the  dispersion  relationship  in  deepwater.  The  value  of 
2(-^^  Is  thus  found  to  be  approximately  equal  to  8.4  x 10~^. 

The  equilibrium  function  in  frequency  domain  can  be  readily  shown  as 


^a^”^  * c"  <36) 

g I k-f^(n) 

where  Cg  is  the  group  velocity  and  k - fj^(n)  is  the  wave  dispersion  relation- 
ship in  shallow  water  • Since,  in  general,  k is  a function  of  both  water 
depth  and  current,  the  equilibrium  function,  or  g^(k)  should  also  vary 

with  these  parameters.  By  virtue  of  Eq.  (10)  which  is  commonly  accepted 
dispersion  relationship,  the  variation  of  ^^Cn)  or  g^(k)  with  respect  to  both 
u and  d can  be  easily  established.  Figures  32  and  33  plot  the  ^^(n)  as  a 
function  u keeping  d constant.  Figures  34  and  35,  on  the  other  hand,  shows 
the  variation  of  4>^(n)  with  respect  to  d for  a constant  u.  It  can  be  seen 
that  both  u and  d (for  u 0),  in  particular,  u,  affect  slenif icantly  the 
equilibrium  function.  A counter  current,  for  instance,  will  surpass  a signifi- 
cant portion  of  energy  in  the  high  frequency  end  and  there  is  a cut  off 
frequency  beyond  which  the  energy  is  zero.  For  the  case  of  a null  current, 

shown  to  be  invariant  with  respect  to  water  depth;  g (n),  however. 


While  the  as  derived  seems  to  explain  reasonably  well  the  behavior 

of  high  frequency  components  in  shoaling  water,  it  also  indicates  that  for  small 
kd  (long  wave  components  in  shallow  water),  waves  can  hardly  reach  saturation 
and  thus  no  breaking  will  occur  in  this  range.  The  latter  result  indeed  is 
puzzling  and  contradicts  to  observation.  This  problem  arises  from  the  basic 
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assumption  in  that  the  breaking  criterion  is  based  upon  vertical  acceleration. 
As  waves  enter  into  shallow  region,  the  obitor  water  particle  mtion 
becomes  elliptic.  Thus,  Che  fluid  becomes  more  and  more  difficult  to  attain 
the  critical  vertical  acceleration  to  reach  incipient  breaking.  While  the 
critical  vertical  acceleration  might  still  remain  as  a sufficient  condition, 
it  no  longer  is  necessary  and  the  wave  might  break  when  other  modes  of 
motion  exceed  certain  critical  value.  For  wave  components  of  small  wave 
numbers  in  shallow  water  the  fluid  motion  in  the  horizontal  direction,  both 
velocity  and  acceleration,  becomes  significantly  larger  than  the  vertical 
one.  Therefore,  waves  may  break  when  the  horizontal  fluid  motion  exceeds 
certain  value  far  before  the  vertical  acceleration  reaches  the  breaking 
criterion.  The  physical  appearances  of  these  two  kinds  of  breaking  might, 
also,  be  quite  different;  the  fluid  separates  from  the  wave  in  the 
horizontal  direction  in  one  case  and  in  the  vertical  direction  in  the 
other.  They  both,  however,  can  serve  as  effective  means  of  restraining 
the  wave  growth. 

A commonly  invoked  criterion  in  the  horizontal  direction  is  that  the 
horizontal  fluid  velocity  at  the  surface  exceeds  a certain  fraction  of  wave 
celerity.  Again  utilizing  linear  wave  theory,  this  criterion  can  be 
expressed  as 


(37) 


where  u^  is  the  horizontal  component  of  water  particle  velocity,  Oj^  is  a 
constant  and  c is  the  wave  celerity.  For  a monochromatic  wave  train,  one 
may  argue  from  physical  intuition  that  a wave  will  break  when  the  horizontal 
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fluid  motion  exceeds  the  wave  celerity,  or  equivalently,  <x^  should  be  equal 
to  unity.  Field  observation  and  experimental  evidences  seem  to  indicate 
that  should  be  around  0.5. 

The  widely  used  Miche*s  (1944)  limiting  steepness  wave  in  shallow 
water,  for  Instance,  can  be  shown  as  equivalent  to  Eq.  (37)  with  a value 
of  “ 0,45.  The  modified  Mlche*s  criterion  as  used  in  this  work  (Eq.  21) 
will  yield  a value  of  * 0.38.  It  is  also  Interesting  to  observe  here  that 
the  criterion  expressed  in  the  above  equation  when  applied  to  deepwater 
becomes  equivalent  to  the  acceleration  criterion  of  Eq.  (33),  provided 
a - a^,  Longuet-Higglns  (1976)  recently  showed  that  the  wave  will  break 
when  attains  0,397  g,  in  other  words,  u “ 0,397  which  is  Indeed  compatible 
with  the  range  of  values  of  indicated  here. 

If  the  same  procedures  used  to  arrive  at  the  acceleration-limited 
equilibrium  function  are  followed,  the  velocity-limited  counterpart  can  be 
obtained.  It  assumes  the  following  form: 


f 

I 


f 


“l^2  tanh  kd  (tanh  kd  - l«l  sech  kd) 


When  matched  with  the  deepwater  equilibrium  condition,  the  coefficient 

-2 

is  found  to  be  equal  to  9.4  x 10  . 


(38) 


Like  G (k). 


[ 

t 


normalized  function  G^(k)  can  be  defined  such  that 


G^(k) 


(39) 
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Similar  to  the  acceleratlon-^llmlted  equilibrium  function^  the  value  of  k^g  (k) 

3 

or  k G^(k),  Is  a function  of  kd  only*  This  relationship  is  plotted  In  the 

same  graph  with  k^  ^^8-  30  to  offer  comparison.  In  contrast  to 

3 3 

k G^(k),  the  function  k G^(k)  is  a monotonlcally  Increasing  function  with 

3 

increasing  kd  and  approaches  asymptotically  to  the  constant  value  as  k G^(k) 
when  kd  -►  *, 

3 

Field  data  seemed  to  support  the  trend  as  depicted  by  the  curve  k G^(k). 
Pierson’s  data  from  the  stereo'wave  observation  project  (1962)  is  plotted  in 
the  same  figure  to  indicate  this  trend*  In  plotting  these  data,  an  arbitrary 
water  depth  d " 10  m was  assumed.  Thus,  the  comparison  is  only  qualitative. 

The  variations  of  G^(k)  with  kd  is  shown  in  Fig,  31  to  offer  a comparison  with 
Ga(k).  Unlike  G^Ck),  G^(k)  assumes  a bell-shape  and  has  a maximum  at  kd  = 0.072. 
This,  at  least,  offers  a possible  explanation  as  to  why  most  of  the  energy 
spectra  obtained  in  the  field  are  bell^shaped  and  have  a selected  range  of  wave 
number  or  frequency  that  has  a larger  energy-containing  capacity. 

The  corresponding  equilibrium  function  in  frequency  domain  based  upon 
velocity  criterion,  if^Cn)  can  be  computed  in  a similar  fashion  as  ♦g(n).  Unlike 
$^(n),  is  not  an  invariani  with  d when  u * 0.  Examples  of  ♦^^(n)  as 

function  of  u with  d * constant  are  plotted  in  Figs.  32  and  33  with  the 
corresponding  ♦^(n)  curves.  The  variation  of  with  d keeping  u 

constant  is  shown  in  Figs.  34  and  35  where  <^(n)  curves  are  plotted  with  the 
same  condition.  From  these  figures,  it  can  be  seen  that  <^^(n)  and 
intersect  at  a certain  frequency  n^.  For  wave  components  of  n > n^,  the 
equilibrium  is  maintained  by  acceleration  criterion.  When  fluid  acceleration 
exceeds  this  criterion,  breaking  occurs  with  the  appearance  of  spilling 


3 

i 


- 257  - 

over*  For  wave  components  of  n < n^,  the  equilibrium  is  governed  by  the 
horizontal  fluid  velocity.  When  fluid  velocity  exceeds  this  limit, 
breaking  also  occurs  with  the  appearance  of  curling  over,  or  commonly  known 
as  plunging.  This  composite  picture  is  Illustrated  in  Fig.  36. 


A number  of  comparisons  between  field  data  and  the  composite  equilibrium 
curves  are  shown  In  Figs.  37  to  38  . The  cases  selected  here  are  either  fully 
breaking  or  partially  breaking  waves.  The  current  used  In  these  cases  are  not 
actually  measured  values  but  arbitrary  assigned  to  yield  better  comparisons. 
The  range  of  these  values  are,  however,  not  unreasonable  under  the  actual 
conditions. 

8.  Wave  Climate  at  the  Test  Site 

Based  on  field  observation,  aided  by  the  computer  prediction,  some 
assessment  can  be  made  on  the  wave  climate  at  the  test  site.  Since  the  field 
measurement  was  conducted  for  a relatively  short  period  with  the  highest 
measured  significant  wave  height  up  to  only  1.5  M neither  statistical  Infor- 
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nation  can  be  derived  nor  the  results  can  be  extrapolated  with  any  confidence 
to  waves  of  much  larger  amplitude.  In  addition,  the  measurement  was  obtained 
on  a single  line  and  the  prediction  covered  a section  approximately  400  m from 
each  side  of  the  control  line;  the  observations  made  herewith  should  also  be 
viewed  within  this  context. 

The  shoreline  at  the  test  section  forms  an  angle  of  approximately  17"  from 
the  north  as  shown  In  Fig.  7.  The  offshore  bar  serves  as  an  effective  wave 
breaker  during  low  tide.  The  breaker  role  apparently  Is  Insignificant  during 
high  tide  for  waves  of  common  occurrence.  However,  even  at  high  tide,  the 
bar  seems  to  act  as  a wave  reflector  and  offers  partial  protection  when  waves 
are  approaching  from  S-W  quadrant.  Waves  from  N-W  quadrant  are  manifestly  less 
affected  by  the  presence  of  the  bar.  If  no  breaking  occurs  at  the  bar,  waves 
from  Che  N-W  quadrant  are.  In  general,  amplified  when  approaching  shoreline; 
the  wave  amplification  as  well  as  the  Incident  angle  appears  to  be  rather  uniform 
along  the  test  section;  the  energy  density  function  peaks  up  In  the  vicinity 
of  dominant  frequency  components.  This  situation  Is  Illustrated  In  Table  3 for 
the  wave  component,  n - 0.15.  For  waves  approaching  from  S-W  quadrant,  the 
situation  Is  quite  different.  The  variations  of  both  wave  angle  and  wave  height 
are  pronounced  from  grid  to  grid  In  both  longshore  and  offshore  directions.  The 
energy  density  function,  although  peaks  up  In  general  when  close  to  shoreline 
may  actually  diminish  at  certain  localities.  This  situation  Is  Illustrated  In 
Table  4 for  the  wave  component  n - 0.15.  Fig.  39  compares  the  spatial  variations 
of  near-breaking  wave  angles  for  the  case  when  waves  approach  from  N-W  quadrant 
to  that  from  S-W  quadrant.  It  is  quite  apparent  that  In  the  former  case,  the 
breaking  angles  are  rather  orderly  oriented  whereas  In  the  latter  case,  the 
breaking  angles  are  more  erratic.  Since  wave-induced  currents,  such  as 
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Fig«  39:  Comparison  of  spatial  variations  of  near~breaking  wave  angles 

approaching  from  N-W  quadrant  to  that  from  S-W  quadrant 


the  longshore  drift  Is  governed  by  wave  angle  and  wave  height i one  expects 
that  along  the  test  section  waves  from  the  N-W  quadrant  should  be  owre 
effective  In  creating  an  undlsrupted  longshore  current  than  that  from  the 
S-W  quadrant.  The  latter  offers  a better  chance  to  create  circulation 
cells. 

As  pointed  out,  already  In  previous  sections,  waves  from  the  S-W  are 
more  varied  spatially  and  are  more  sensitive  to  the  bottom  changes.  This 
point  Is  well  Illustrated  In  Fig.  29  where  the  predicted  spectrum  at  75  m 
(the  adjacent  grid)  south  of  station  2 Is  compared  with  that  at  Station  2; 
the  change  Is  quite  significant  for  two  adjacent  grids. 

Although  the  offshore  bar  plays  a protection  role  most  of  the  time, 

It  also  could  occasionally  become  an  energy  trapper  to  prevent  wave  energy 
from  escaping  the  Inshore  zone.  This  occurs  mainly  under  the  offshore 
wind  condition  as  sho%m  in  Fig.  13.  Fortunately,  this  condition  often 
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is  associated  with  moderate  waves,  thus,  may  not  pose  a serious  condition 
to  shore  erosion.  However,  the  energy  piling  up  on  the  shoreward  side  of 
the  bar  might  be  just  sufficient  to  maintain  the  permanence  of  the  bar. 

9.  Conclusions 

The  field  data  obtained  at  the  Island  of  Sylc  supports  the  validity 
of  the  numerical  model  on  wave  energy  transformation  in  regions  with  as 
complex  a nearshore  hydrograph  as  the  test  site.  With  due  consideration 
on  the  complexity  of  the  problem  and  the  simplified  assumptions  made  in 
the  modeling,  the  numerical  predictions  yield  quite  acceptable  results 
of  energy-containing  wave  components  and  long-wave  components.  In  these 
ranges,  the  agreement  between  predicted  results  and  field  data  is  considered 
good  up  and  prior  to  breaking.  It  suggests  that  the  transformation  of 
waves  components  within  these  ranges  Is  mainly  Influenced  by  shoaling 
and  refraction,  the  two  factors  Included  in  the  modeling.  The  numerical 
results  are  unsatisfying  in  high  frequency  components  where  local  effects 
of  energy  dissipation,  generation  or  transfer  are  deemed  to  be  more 
Important.  All  these  effects  are  not  being  considered  in  the  present 
model. 

Turbulent  generation  due  to  wave  instability  seems  to  be  the  dominant 
mechanism  of  energy  dissipation  as  opposed  to  bottom  friction  which  seems 
to  play  a very  minor  role  except  in  the  region  of  wave  uprush.  Based  upon 
commonly  used  wave  Instability  criteria,  the  equilibrium  energy  spectral 
density  function  in  shoaling  water  has  been  developed.  This  function 
provides  saturation  condition  on  spectral  density  for  the  complete  frequency 
(or  wave-number)  range  and,  thus,  enables  one  to  estimate  the  limiting 


- 264  - 

spectra  as  functions  of  local  water  depth  and  current  conditions. 

Comparisons  with  field  data  are  encouraging,  in  particular,  for  cases  where 
waves  break  over  offshore  obstacles  such  as  the  offshore  bar  in  the 
present  study.  An  explanation  has  been  offered  on  the  difference  between 
spilling  and  plunging  breakers. 

The  limited  field  data  alone  does  not  provide  sufficient  Inforaatlon 
to  perform  wave  climate  assessment  at  the  test  location.  However,  when 
coupled  with  numerical  computations  which  extend  the  point  infonnat^lon  to 
spatial  distribution  certain  general  observations  can  be  made.  The  offshore 
bar  plays  the  double  roles  as  energy  dlsslpator  and  trapper  depending  mainly 
on  the  wind  direction, tidal  stage  and  wave  height.  For  low  tide,  the  role  of 
dissipation  dominates.  For  high  tide  and  nonbreaking  wave,  the  bar  acts 
like  a partial  reflector  under  onshore  wind  and  an  energy  trapper  under  off- 
shore wind.  Waves  fron  the  N-W  quadrant  are  found  to  be  less  affected  by 
the  bar  and  are  more  effective  in  generating  unidirectional  longshore  drift 
directing  southward.  Waves  from  the  S-W  quadrant  appear  to  have  greater  spatial 
variations  both  in  wave  angle  and  wave  magnitude  and  are  likely  to  promote  nearshore 
circulation  cell.  For  wind  wave  conditions,  the  JONSWAP  spectra  are  reasonable 
deepwater  Inputs  for  the  tested  region.  All  the  above  observations  await  further 
field  experiment  for  verification. 

In  conclusion,  the  numerical  model  has  shown  to  be  valid  and  provides  a 
sueful  tool  to  yield  spatial  information  on  wave  energy  spectra  In  nearshore 
zone.  Further  improvements  could  be  made  in  the  following  directions: 

1.  Field  data  with  simultaneous  current  information. 

2.  The  effect  of  directional  spreading. 

3.  The  Interaction  among  wave  components 


I 
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and  in  that  order  of  importance*  Analytical  and  experimental  work  should 
be  continued  in  Che  area  of  developing  limiting  wave  spectra  in  shallow 
water.  This  information  should  be  useful  in  determining  limiting  design 
condition  for  engineering  purposes. 
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Appendix  1 

Computer  Progreun-Periodogram  Generator 
and  Spectral  Analysis 


rKrmfw«. 


275 


o o 


o 


^ *n  i/i  c/>  O 

C «««««« 

• Ui  MW  ww  ww 

<«  9 s 9-^  o<«  ca« 

■ aa  • K.  1^ 

•»*»-0©*^«-‘NJl>JPO»n'03l/>»»«-0»/»rNjrMOl/»*0»no</»3 
♦ ♦WM  •WWMWMM  Z»^WM««MM««MM««Z 
»^^ww^wwwwwwO*«  w«k>WMwwWMw%,»MWm 

■ Arf  K>  O p*fc<VAI^^F>^^0 


irt 


>a 


K-  o W fcO 
ui  a I I 

C (A  c UJ  < 

I W a. 

Z « O < 

uj  z X : 

^ M (/)  I 
Z Ui  U.  I 

O C V Ik  I 
^ HO 
H Z « 
w V)  H • 
Ui  Z < 
Z H Ui 
UI  cd  o X • 
ct  oe  H : 
Ui  ^ X I 
>•  Z < UJ  I 
Ui  X o I 
Z H 
• K W X I 


O H 
V)  z Z 

Z U.  Ui 

Oc  H 

Ui  z (/» 


Z Ik  < 

O O O Z 3 

■HO  o : 

I z H M : 

Z o « H 

t < H OS  Z J 


I o o ec 

■ Z H </)  Ui  I 

I z H U4  H I 


-I  O H 

> z « z o 

• u«  Z Z 

: z o M z 

> VO  H M 


» O fi 

z ^ 

Ui  ^ H 
H >0  «/> 

Z S3 

e z 

H 

»>  O Ui 
Vi  ••  >o  -J 
Z ••  s VO 
Ui  — X z 
3 ..  Kt  < 

® C 
Ui  in  z 
Z </>  H 

Ik  Ui  i/>  • 

H H Z X 
Ui  O < H 


I Z O Z < * 
V sm  « H 
> IM  S < < « 

I M o so: 

' Z lA  A 

-3  H : 
M >c  w>  3 
Z o Z i 

S S 3 UI  Z 

. VO  A Z i-u  i 


' Z Z (B  V 
• O O Z 

I r X Ui  z 

I X z z < 

I O O M Ui 

sw  w z or 


Ui<.iCiWWWMViUiOCiWU*Wi.^CiV.>V,iWt.iU>OCl 


z o 

O A 


^ Z 


II  «• 

-»  H 

% M 


« X 
S ►“ 
a 

o (/> 


« (/)  X O X 
K a O < h> 


^ o 

o ^ 

o 

O 'O 


O C9 
<NJ  K» 
3 O 
'C 


U>  K Z % 

o o « %» - 


^ ^ 


) «o  ^ 

« « II 

> >0  m 


O O 
>0  K 

o o 

'O  <0 


O O </)  N.  M 


wuiujuJUiUJOO^Uiw 


: a « oe  de  o 


-*  II 
O -1 

M 

O 

a a 


oc  oe  oe  or  O 


fvj  o a. 

^ • •<«.  o o 

^ 1/^  rg  9 • • 

iM^O  • rg  fvj  r*  O 

'O  •«  c « « o < 

00  •o.v^aanvs 
• Xb  M II  II  II  II  0 
II  rg  a rsi  r>j  (/>  u 
II  MM«M>>fiLCO 
tBO.a.oraauia 


or  X 
^ ^ 
tt  U H> 
a uj  9 

*/»  a.  o 

* tA  at 
» CD 
AJ  Ui  S 
CO  >-  t/» 

• 3 

<M  O.  (/) 
H C M 

a o X 


w o 

Z M 
X 

I-  ^ ♦ 

CL  II  «*. 

uj  X a 


^ A 

« o 

u»  a (A 


M ffi  I 

a.  3 : 


<N  M 
II  X 


I O 4U 
■ -J  at  ^ 
I Ui  s A 
a 3 Z 3 z 
O (/)  tA  </>  < 

_l  (/)  _i  tA  _J 
^ ^ M _J 


* 

c s « « 

% ^ ^ ^ 

•-  »“  II  II 

II  II  « 

^ ^ * z 

Z Z \ « 

^ ^ II  II 

»l  II  -»  ^ 

“>-»%% 


300000  OOC  (OOOO 
r*N.K*N.'#iy^Ofv(M 

'O'A'O'O'C'O'O'O'O  ▼*<C'0'0 

I'O'O'O'O'O'O'O’O'OOM'^'O'OO 

wwwz 

«uJWujUiuJiuuiuiUiHi40uJuJUi*H 


3 0 3 0 0 3 
v ry  r<^  ^ rsj 

O C O O O O O 
l/>  UN  t/N  •/)  'C 


AD-A039  988  DELAWARE  UNIV  NEWARK  DEPT  OF  CIVIL  ENGINEERING  F/G  8/3 

MEASUREMENTS  AND  COMPUTATION  OF  WAVE  SPECTRAL  TRANSFORMATION  AT—ETC(U) 
NOV  76  H WANGr  W YANG  N00014-76-C-0392 

OCEAN  ENGINEERING- 11  NL 


UNCLASSIFIED 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  Of  STANDARDS  1963  A 


290 


1 


CONTINUE 

IF  (1FLA6.E0.1  ) &0  TO  10 


\ 


i 


1 


ZusanuLenfassurg 

-iessungen  und  Berechnungen  der  Veranderur4g 
von  Wellenspektren  vor  der  Insel  Sjlt 

Im  Kiistenvorfeld  der  Insel  Sylt,  rd.  90C  m seewarta,  wurden  Wellen- 
spektren  geoessen,  Die  im  Mai  1576  fUr  die  Dauer  von  zwei  Tocher 
betriebene  iieBkette  urnfaBte  sieben  Stationen*  Die  Srgebnisse  der 
Naturmesaungen  wurden  nach  der  Auswertur.g  oit  einem  von  den  Ver- 
fassern  fiir  dieses  Projekt  vorher  entwickelten  numerischen  Berech- 
nungsansatz  fur  die  2nergieumwandlung  in  Flachwassergebieten  ver- 
glichen* 

Die  Terte  der  numerischer  Berechnung  zeigten  eine  gute  t*berein- 
stifflinung  ait  den  Daten  aua  der  Natur  is  Bereich  der  unteren  und 
mittleren  Frequenzeni  im  Bereich  hoher  Frequenzen  war  die  tberein- 
stioffiung  weniger  gut* 

Ausgehend  von  den  bekannten  Stabilitatsk  ‘iterien  fur  Tellen  wurde 
eine  Dleichgowichtsfunktion  dvr  spcktralen  Lnergiedichte  fiir  Berei- 
che  mit  beachrankter  Tassertiefe  entwlckelt,  die  fiir  das  gesamte 
Spektrum  gute  Vhereinstimmung  hinsichtlich  der  Spektraldichte  zeigt. 

Die  Vergleiche  mit  den  Naturoessungen  fallen  ebenfalls  gut  aus,  es 
bedarf  jedoch  weiterer  kessungen  fur  eine  allgemeingul tige  3e- 
statigungo 

Aus  der  Verbindur.g  von  MeSergebnissen  und  theoretischen  Berechnun- 
gen  konnten  Aussagen  iiber  das  V/ellenklioa  des  Untersuchungsgebie- 
tes  vor  Sylt  gemacht  werden.  Es  wurde  gezeigt,  daC  das  kUstenparal- 
lele  Riff  eine  Doppelfunktion  als  linergie  utwandelnde  und  Energie 
konzentrierende  uorphologische  Formation  hat*  Teiterhin  konr.te  ab- 
geleitet  werden,  daS  Wellen  aus  dem  ir.V-  'Auadranten  den  groflten  Ein- 
flu£  auf  die  Entstehar.g  von  Kustenlangsstromangen  und  damit  auf  ^ 

den  Xiistenlangstransport  haben.  J 
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